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Abstract. We present analysis for the reduction of an inertial description of fluid-structure in- 
teractions subject to thermal fluctuations. We show how the viscous coupling between the immersed 
structures and the fluid can be simplified in the regime where this coupling becomes increasingly 
strong. Many descriptions in fluid mechanics and in the formulation of computational methods ac- 
count for fluid-structure interactions through viscous drag terms to tranfer momentum from the fluid 
to immersed structures. In the inertial regime, this coupling often introduces a prohibitively small 
time-scale into the temporal dynamics of the fluid-structure system. This is further exacerbated in 
the presence of thermal fluctuations. We discuss here a systematic reduction technique for the full 
inertial equations to obtain a simplified description where this coupling term is eliminated. This 
approach also accounts for the effective stochastic equations for the fluid-structure dynamics. The 
analysis is based on use of the infinitesmal generator of the SPDEs and a singular perturbation anal- 
ysis of the backward kolomogorov PDEs. We also discuss the physical motivations and interpretation 
of the obtained reduced description of the fluid-structure system. 
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1. Introduction. The investigation of fluid-structure interactions has a rich history. 
Problems motivating past work include the study of aerodynamic oscillations induced in 
airplane wings and propellers [20, 22], the study of animal locomotion including swimming 
and insect flight [35, 29, 41], and the study of physiological problems such as blood flow 
through heart valves [25, 39, 23]. Central to studying these problems in practice is the 
development of descriptions which capture essential features of the fluid structure interactions 
while introducing approximations which facilitate analysis and the development of tractable 
numerical methods [20, 36]. This area of research is still very active [25, 11, 36, 10]. While 
we have above mentioned macroscopic problems involving fluid-structure interaction, recent 
developments in the sciences and technology motivate the study of fluid-structure interactions 
in new physical regimes often involving very small length and time scales. This includes the 
dynamic responses of soft materials and complex fluids [44, 46, 13, 37]. At the molecular 
level advances in computational modeling at the level of coarse-grained descriptions with 
implicit solvent also motivate the need for the development of new descriptions of fluid- 
structure interactions to account through continuum fields the important roles played by 
momentum transfer of the neglected solvent fluid in dynamic responses. This includes the 
study of implicit solvent models for polymers, gels, and lipid bilayers (cite). For these 
problems thermal fluctuations often play an important role and pose additional challenges 
in the study of fluid-structure systems. 

Significant past work has been done on the formulation of descriptions of fluid-structure 
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interactions subject to thermal fluctuations. Many of these analytic and numerical ap- 
proaches originate from the polymer physics community [17, 21, 42, 9]. To obtain descrip- 
tions tractable for analysis and numerical simulation, these approaches typically place an 
emphasis on approximations which retain only the structure degrees of freedom. This of- 
ten results in significant simplifications in the descriptions and in significant computational 
savings. This eliminates the many degrees of freedom associated with the fluid and avoids 
having to resolve the potentially intricate and stiff stochastic dynamics of the fluid. These 
approaches have worked especially well for the study of bulk phenomena in free solution and 
the study of complex fluids and soft materials [42, 17, 28]. 

Recent applications arising in the sciences and in technological fields present situations in 
which resolving the dynamics of the fluid may be important and even advantageous both for 
modeling and computation. This includes modeling the spectroscopic responses of biological 
materials [47, 24, 34], studying transport in microfluidic and nanofluidic devices [44, 38], 
and investigating dynamics in biological systems [2, 15]. There are also other motivations 
for representing the fluid explicitly and resolving its stochastic dynamics. This includes the 
development of hybrid fluid-particle models in which thermal fluctuations mediate important 
effects when coupling continuum and particle descriptions [16, 18], the study of hydrodynamic 
coupling and diffusion in the vicinity of surfaces having complicated geometries [44] , and the 
study of systems in which there are many interacting mechanical structures [6, 40, 39, 3, 8, 
10, 45, 36, 48, 26, 31, 5, 12, 7, 19, 6]. To help obtain less stiff equations for amenable to 
practical numerical approximation, some prior work has been done on the above mentioned 
fluid-structure interaction approaches. For example, for the Stochastic Immersed Boundary 
Method a reduction procedure was developed to eliminate the fluid degrees of freedom and 
obtain a purely particle based description [32]. 

To facilitate the development of tractable descriptions and efficient numerical methods 
for study of fluid-structure systems, we consider a rather general formalism which captures 
essential features of the coupled stochastic dynamics of the fluid and structures. To model 
the fluid-structure system, a mechanical description is utilized involving both Eulerian and 
Lagrangian reference frames. Such mixed descriptions arise rather naturally, since it is often 
convenient to describe the structure configurations in a Lagrangian reference frame while it 
is convenient to describe the fluid in an Eulerian reference frame. In practice, this presents 
a number of challenges for analysis and numerical studies. A central issue concerns how to 
couple the descriptions to represent accurately the fluid-structure interactions, while obtain- 
ing a coupled description which can be treated efficiently ultimately by numerical methods. 
An important issue often arising concerns the challenge of developing temporal integrators 
capable of handling a wide range of time-scales exhibited by the fluid-structure dynamics. 
Another important issue concerns how to account properly for thermal fluctuations in such 
numerical integrators and associated approximations. This all must be done carefully to be 
consistent with statistical mechanics. 

To address these challenges and gain further insights into the fluid-structure equations, 
we develop a systematic reduction of the full stochastic description to a simplified description 
that removes the fastest temporal dynamics of the system and reduces the number of degrees 
of freedom. Our approach is based on the Infinitesmal Generator of the SELM SPDEs and 
a singular perturbation analysis of the Backward Kolomogorov PDEs following an approach 
recently introduced in [33, 32]. For the SELM approach, we consider the regime in which 
the fluid-structure coupling becomes strong. We then discuss our reduced descriptions of the 
fluid-structure interactions in a few different regimes and their physical interpretation. 

2. Summary of the Stochastic Eulerian Lagrangian Method. We sum- 
marize here the Stochastic Eulerian Lagrangian Method, abbreviated as SELM. We present 
the general formalism and a number of alternative descriptions of the fluid-structure system. 
In many situations the stochastic differential equations for the full fluid-structure dynamics 
exhibits stiffness. To cope with this issue and to develop efficient numerical methods, sim- 
plified descriptions are discussed for various physical regimes. A more detailed discussion 
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Fig. 2.1. T/ie description of the fluid- structure system utilizes both Eulerian and Lagrangian 
reference frames. The structure mechanics are often most naturally described using a Lagrangian 
reference frame. The fluid mechanics are often most naturally described using an Eulerian reference 
frame. The mapping X(q) relates the Lagrangian reference frame to the Eulerian reference frame. 
The operator T prescribes how structures are to be coupled to the fluid. The operator A prescribes 
how the fluid is to be coupled to the structures. A variety of fluid-structure interactions can be 
represented in this way. This includes rigid and deformable bodies, membrane structures, polymeric 
structures, or point particles. 



and derivation of SELM and the reduced equations in each of the physical regimes is given 
in Section 6. 

To study the dynamics of fluid-structure interactions in the presence of thermal fluc- 
tuations, we utilize a mechanical description involving Eulerian and Lagrangian reference 
frames. Such mixed descriptions arise rather naturally, since it is often convenient to de- 
scribe the structure configurations in a Lagrangian reference frame while it is convenient to 
describe the fluid in an Eulerian reference frame. In principle more general descriptions us- 
ing other reference frames could also be considered. Descriptions for fluid-structure systems 
having these features can be described rather generally by the following dynamic equations 

=£u + A[T(v-ru)]+A + f thnl (2.1) 

dv 

m— = -T (v - Tu) - Vx$[X] + C + F thm (2.2) 



dt 
rfX 
~dt 



(2.3) 



The u denotes the velocity of the fluid, p the uniform fluid density. The X denotes the 
configuration of the structure and v the velocity of the structure. The mass of the structure 
is denoted by m. To simplify the presentation we treat here only the case when p and m 
are constant, but with some modifications these could also be treated as variable. The A, £ 
are Lagrange multipliers for imposed constraints, such as incompressibility of the fluid or a 
rigid body constraint of a structure. The operator C is used to account for dissipation in the 
fluid, such as associated with Newtonian fluid stresses [1]. To account for how the fluid and 
structures are coupled, a few general operators are introduced, V, T, A. 

The linear operators V, A, T are used to model the fluid-structure coupling. The V 
operator describes how a structure depends on the fluid flow while — T is a negative definite 
dissipative operator describing the viscous interactions coupling the structure to the fluid. 
We assume throughout that this dissipative operator is symmetric, T = T T . The linear 
operator A is used to attribute a spatial location for the viscous interactions between the 
structure and fluid. The linear operators are assumed to have dependence only on the 
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configuration degrees of freedom T = T[X], A — A[X]. We assume further that T does not 
have any dependence on X. 

To account for the mechanics of structures, $[X] denotes the potential energy of the 
configuration X. The total energy associated with this fluid-structure system is given by 

B[u,v,X] = J ip|u(y)| 2 dy + l -mv 2 + $[X]. (2.4) 

The first two terms give the kinetic energy of the fluid and structures. The last term gives 
the potential energy of the structures. 

As we shall discuss, it is natural to consider coupling operators A and Y which are 
adjoint in the sense 

/ (ru)(q)-v(q)dq= / u(x) • (Av)(x)dx (2.5) 
J s Jn 

for any u and v. The S and denote the spaces used to parameterize respectively the 
structures and the fluid. We denote such an adjoint by A = r 1 " or F — A*. This adjoint 
condition can be shown to have the important consequence that the fluid-structure coupling 
conserves energy when T — s> oo in the inviscid and zero temperature limit. 

In practice, the conditions discussed above can be relaxed somewhat. For our present 
purposes these conditions help simplify the presentation. Each of these operators will be 
discussed in more detail below. 

To account for thermal fluctuations, a random force density f thm is introduced in the 
fluid equations and F t h m in the structure equations. These account for spontaneous changes 
in the system momentum which occurs as a result of the influence of unresolved microscopic 
degrees of freedom and unresolved events occurring in the fluid and in the fluid-structure 
interactions. 

The thermal fluctuations consistent with the form of the total energy and relaxation 
dynamics of the system are taken into account by the introduction of stochastic driving 
fields in the momentum equations of the fluid and structures. The stochastic driving fields 
are taken to be Gaussian processes with mean zero and with 5-correlation in time [43]. By 
the fluctuation-dissipation principle [43] these have covariances given by 

(tthm(«)C»(*)) = - ( 2k BT) (C - ATr) 5(t - s) (2.6) 
(Fthm(*)F£»(t)) = (2k B T) TS(t - s) (2.7) 
(fthm(*)FS»(t)) =-(2* B T)AT*(*-«). (2.8) 

We have used that T — A 1 " and T = T T . We remark that the notation gh T which is used 
for the covariance operators should be interpreted as the tensor product. This notation is 
meant to suggest the analogue to the outer-product operation which holds in the discrete 
setting [4] . A more detailed discussion and derivation of the thermal fluctuations is given in 
Section 6. 

It is important to mention that some care must be taken when using the above formalism 
in practice and when choosing operators. An important issue concerns the treatment of the 
material derivative of the fluid, du/dt — du/dt + u • Vu. For stochastic systems the field 
u is often highly irregular and not defined in a point-wise sense, but rather only in the 
sense of a generalized function (distribution) [14, 30]. This presents issues in how to define 
the non-linear term arising in the material derivative, which appears to require point-wise 
values of u. For such irregular velocity fields, this also calls into question the applicability of 
the theorems typically used to derive the differential equations from the conservation laws. 
For instance, for such velocity fields the fluid material body may no longer exhibit smooth 
deformations over time. 

There are a number of ways to deal with this issue. The first is to consider a regular- 
ization of the fluid stresses, which are typically the source of irregularity, see equation 2.6. 
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This can be motivated by the fact that the fluid stress tensors typically considered in con- 
tinuum mechanics are expected to become inaccurate at molecular length-scales. Ideally, 
from molecular models of the fluid the small-length scale (large wave-number) responses of 
the fluid could be determined and provide a justified regularization. For instance, this could 
provide an alternative to using responses based on Newtonian stresses for all length-scales. 
For the SELM formalism, this would simply correspond to using for C an alternative to 
the dissipative operator based on Newtonian stresses. The second more easily implemented 
approach is simply to work with the linearized material derivative, which still retains many 
of the essential features of the fluid dynamics and is useful for many applications [5] . 

In this version of SELM, we assume the advective derivative is well-defined and include 
its terms in the computation. We will use a time scale specific to the advective term to keep 
track of which terms arise due to the convective term. In this way, the advective term may 
be neglected in specific regimes where the approximation du/dt — du/dt can be used. 

3. Formulation in terms of Total Momentum and Particle Momentum 
equations. We now consider the regime in which the full dynamics of the fluid-structure 
system are retained but reformulated in terms of a field describing the total momentum of the 
fluid-structure system at a given spatial location. This description is more convenient to work 
with in practice since it results in simplifications in the stochastic driving fields. Primarily, 
the two noise terms in the stochastic differential equations involved become uncorrelated, 
simplifying the analysis. For this purpose we define 



p(x(<),t)=/ni(x(t),t) + A(X(t))[mv(<)]. (3.1) 

The operator A is used to give the distribution in space of the momentum associated with 
the structures for given configuration X(t). We may take a time derivative in p to obtain: 

dp du „ , . . . dv . ,„ „. 

Now, 

= + P u ■ Vu = £u + A[T(v - m)] + A + f thm (3.3) 
Using this approach, the fluid-structure dynamics are described by 

^=£u + A[T(v-ru)]+A + f tta (3.4) 

- pu • Vu + V x A(mv) • v + A(m— ) (3.5) 

= £p _1 [p - A(mv)] - p _1 [p - A(mv)] ■ V[p - A(mv)] + V x A(mv) ■ v (3.6) 
+ A[-V X *(X)] + A + Sthm (3.7) 

' /V = -T (v - p" 1 r(p - Amv)) - V X *[X] + C + F thm (3.8) 



dt 
dX 

—r~ = v 
dt 



(3.9) 



where g t i lm = fthm + A[F t hm]- The third term in the first equation arises from the 
dependence of A on the configuration of the structures, A[mv] = (A[X])[mv). The Lagrange 
multipliers for imposed constraints are denoted by A, For the constraints, we use rather 
liberally the notation with the Lagrange multipliers denoted here not necessarily assumed 
to be equal to the previous definition. The stochastic driving fields are again Gaussian with 
mean zero and (^-correlation in time [43]. The stochastic driving fields have the covariance 
structure given by 

<g thm (A»g t T lm (V^)> = - (2k B T)L5 a (X - Y)S(t - s) (3.10) 
(F thm (s)FL 1 (i)) = (2k B T) T S(t - s) (3.11) 
(g*»(X,s)FL(t)) = 0. (3.12) 
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This formulation has the convenient feature that the stochastic driving fields become inde- 
pendent. This is a consequence of using the field for the total momentum for which the 
dissipative exchange of momentum between the fluid and structure no longer arises. In the 
equations for the total momentum, the only source of dissipation remaining occurs from the 
stresses of the fluid. This approach simplifies the effort required to generate numerically the 
stochastic driving fields and will be used throughout. 

4. Non-Dimensional Equations. Our interest in this section is to develop an 
appropriate dimension analysis to make precise the reduce of our stochastic equations in a 
limiting regime corresponding to taking the viscous coupling coefficient large, " Y — > oo" . Our 
reduction procedure will then aim to eliminate equation 3.8, which carries T, and replaces 
equations 3.4 and 3.9 with effective equations in this limiting regime. 

4.1. Dimension Analysis of the Equations. To faciliate proper handling of the 
reduction of the SELM SPDEs, it is convenient to work with non-dimensional equations. 
This will have the added benefit when taking limits of precisely defining the asymptotic 
regime. Since the SELM SPDEs describe physical laws the behaviors of the underlying dy- 
namical system must be invariant under any change of the fundamental units, in this case the 
units of mass, length, time, and temperature. The Buckingham II theorem provides a useful 
guideline for the number of non-dimensional groups required to characterise distinctly the 
non-dimensionally equivalent dynamical systems, in particular, that for TV parameters and 
m fundamental units there is a need for M — N — m non-dimensional constants IIi , . . . IIm 
(cite). Simple counting of the parameters for the SELM SPDEs gives N = 6 and the number 
of fundamental units involved are m = 4 requiring M = 2 non-dimensional groups to be 
identified. 

Before we identify our non-dimensional groups, we will identify the time-scales that will 
allow us to understand the behavior of the system. We will first examine T in 3.8. Suppose 
we can non-dimensionalize T = YoY, where To is a scalar that carries the units of T and 
T is a non-dimensional operator. We will define t\ as the relaxation time for v under the 
T operator, which is evidently m/To. This choice of t\ also characterizes the momentum 
transfer rate between the structure and the fluid, since T controls the dissipative force on 
the structure and fluid. We will take fa = y m£ 2 / ksT , which describes the time-scale for 
the particle to diffuse over a distance I. 

For the purposes of our analysis, the non-dimensional groups we we use are k 
and e = ii/fa. In the notation, p is the fluid density, m the excess particle mass over the 
same volume of fluid, and £ the particle radius. In other words, k describes the density ratio 
of the particle and fluid, while e encapsulates the ratio of the two time-scales involved. We 
will fix k for our reduction. 

Now, we assume the dynamics of the system in which we are interested function on the 
time-scale t%, and that t\ « ti- In this way, the fluctuations due to the interaction between 
the fluid and the structure occur on a much faster time scale than the remaining terms in 
the system. 

Since we are only interested in taking e — > 0, we may use a similar scaling for each of 
our variables other than T. We non-dimensionalize each of the variables using £, m, and ti: 



i = * 2 t,X = £X,v = « v = iv,p = p p= ^P,A = A A= ^A,L = C C= jS^C,® = 



2 



We also non-dimensionalize the noise terms, 



g thm (X,s) = g g th . 



F thm (s) = F F th . 



.m — 



.??). 




(4.1) 



(4.2) 
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Here, 



{ttX,s)^(Y 1 i))=6 3 (X-Y)6 3 (s-i), (4.3) 

(vm T (t)) =5 3 (s~i), (4.4) 

2 ksT . . 

ffo = -^-£0, (4.5) 



2 

Fo 2 = ^T , (4.6) 

12 

2D A D A T = A (4.7) 

where A is a symmetric positive-definite operator. Here, Da is an operator defined with the 
same domain and range as A. Its existence is guranteed by the properties of A. 

4.2. Other Conventions. Since fj(-) and £(X, •) are Gaussian processes with mean 
and ^-correlation, we will be able to use 

V(t)= d §, f (X, t) = *(X) dB ^- , (4.8) 

where 

<0(X)0 T (Y))=<5 3 (X-Y). (4.9) 

Above, the superscript distinguishes the different Weiner processes at each spatial point. 

For the purposes of our analysis, we will approximate p as a function of a finite number 
of spatial points. However, since the number of points we use is arbitrary, the dynamics of 
the system may be approximated arbitrarily well, assuming the system is sufficiently well- 
behaved. Moreover, during our analysis we will use the convention that p = p(t), where we 
establish an implicit spatial dependence of p by using a greater dimension for the range of 
p. For example, if we would like to use N points, we could use (£) where k = 1, N 
corresponds to the spatial point and i is the corresponding component of p at each point. 

With the new convention, there is no need to include cf> because of our discretization 
scheme. Thus the noise due to gth m can be labeled as D-c^j^- 

4.3. Non-Dimensional Equations. Writing each of the terms in their non-dimensional 
form, removing £ and £ for simplicity, and removing the bars from the non-dimensional vari- 
ables for convenience, we can re-write equations 3.4 3.8, and 3.9 as 



dp 

dt 



-KrV[p - A(v)] ■ V x [p - A(v)] + k~ x C (p - A[v]) (4.10) 
A[-Vx*(X)] + (VxA[v]) ■ v + D.c^ (4.11) 



\ [ - Tv + K - 1 Tr(p - A[v]) - eVx$(X)] + ^D r ^ (4.12) 
v (4.13) 



It will be useful to write 4.12 as 



dv 1 „ _ . . / 1 _, dBt , , , , „ 

-=--rC(v-^) + J- e Dr-^. (4.14) 



Here, 



c = (i + « _1 rA), v = c 1 (/t-'rp - eT" 1 v x $(x)) 



(4.15) 
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5. Summary of Results. The result to first order in e is given in 7.78 and 7.79. 
We summarize the results in the limit e — > next. 

The zeroth-order approximation in e is 

^ = -P~Mp ■ V x p - A[mv ] ■ V x p - p ■ V x (A[mv ]) + (A[mv ]) • V x (A[mv ]) (5.1) 
+ mfc B T(V x AC" 1 ) : A T } 

+ p~ x C (p - A[tov ]) + fcsTVxA : C' 1 + (VxA[mv ]) ■ v - AV X $(X) + A + g thm 

f=- < 5 " 2 > 

5.1. Limit with Neglibigle Excess Mass m — > 0. With m — > 0, we note C — > /, 
so that VxA : C^ 1 — » VxA : / = VxA. We obtain the equations 

^ = -p^p ■ V xP ] + p-'Cp + fcsTVx ■ A - AV X $(X) + g thm (5.3) 
f=v (5.4) 

5.2. Regime IV. The description of the fluid-structure system can be further sim- 
plified by considering for the fluid the viscous limit in which p — > oo. In this regime only the 
structure dynamics remain and can be shown to be given by 

^ = # SB lm[-Vx$(X)] + (V x ■ H SELM )k B T + h thm (5.5) 

Hselm = r(-p£)" 1 A (5.6) 
(h tbm («)hS m (t)> = (2fc s T) Hselm 5(t - s). (5.7) 

The p denotes a projection operator imposing constraints, such as incompressibility. The 
adjoint property A = 1^ and symmetry of pC yields an operator H S elm which is symmetric. 
A more detailed discussion and derivation of the equations in this regime is given in Section 6. 

5.3. Summary. This gives an overview of the SELM formalism and the associated 
stochastic differential equations. We remark that each of these regimes were motivated by 
a rather specific limit. Non-dimensional analysis of the equations can also be carried out 
and other limits considered to motivate working with such reduced equations. We discuss 
in more detail the derivation of the reduced equations in our regime of interest in Section 7. 
developed for the SELM formalism in Section 6. 

6. Derivations for the Stochastic Eulerian Lagrangian Method. We now 

discuss formal derivations to motivate the stochastic differential equations used in each of 
the physical regimes. For this purpose, we do not present the most general derivation of the 
equations. For brevity, we make simplifying assumptions when convenient. 

In the initial formulation of SELM, the fluid-structure system is described by 

P^ = Cu + A[T(v - ru)] + A + f thm (6.1) 
= -T (v - Tu) - V X *(X) + C + F thm (6.2) 



dX 

dt 



(6.3) 



The notation and operators appearing in these equations has been discussed in detail in 
Section 2. For these equations, we focus primarily on the motivation for the stochastic 
driving fields used for the fluid-structure system. 
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For the thermal fluctuations of the system, we assume Gaussian random fields with mean 
zero and 5-correlated in time. For such stochastic fields, the central challenge is to determine 
an appropriate covariance structure. For this purpose, we use the fluctuation-dissipation 
principle of statistical mechanics [43, 27]. For linear stochastic differential equations of the 
form 



dZ t = LZtdt + QdBt 
the fluctuation-dissipation principle can be expressed as 



G = QQ 1 



-(LC) - {LCf 



(6.4) 
(6.5) 



This relates the equilibrium covariance structure C of the system to the covariance structure 
G of the stochastic driving field. The operator L accounts for the dissipative dynamics of the 
system. For the equations 6.1 - 6.3, the dissipative operators only appear in the momentum 
equations. This can be shown to have the consequence that there is no thermal forcing 
in the equation for X(i), this will also be confirmed in subsequent sections. To simplify 
the presentation, we do not represent explicitly the stochastic dynamics of the structure 
configuration X. 

For the fluid-structure system it is convenient to work with the stochastic driving fields 
by defining 



q= [p ^tta.m 1 F thll 



(6.6) 



The field q formally is given by q = QdHt/dt and determined by the covariance structure 
G = QQ T . This covariance structure is determined by the fluctuation-dissipation principle 
expressed in equation 6.5 with 



p- L (£ - ATr) 
m _1 Tr 



p _1 AT 
-m _1 T 



C ■■ 



1 k,B TX 



- x k B TX 



(6.7) 
(6.8) 



The X denotes the identity operator. The covariance C was obtained by considering the 
fluctuations at equilibrium. The covariance C is easily found since the Gibbs-Boltzmann 
distribution is a Gaussian with formal density ^(u, v) = ^- exp [— E/UbT]. The Zo is the 
normalization constant for ty. The energy is given by equation 2.4. For this purpose, we 
need only consider the energy E in the case when $ = 0. This gives the covariance structure 



G = (2k B T) 



-p -2 (£-ATr) -m' 1 p- 1 AT 



-i Tr 



rrT z T 



(6.9) 



To obtain this result we use that Y = A' and T = T . From the definition of q, it is found 
the covariance of the stochastic driving fields of SELM are given by equations 2.6- 2.8. This 
provides a description of the thermal fluctuations in the fluid-structure system. 

7. Stochastic Reduction Procedure. 

7.1. Set-Up and Summary of Derivation. The next step in our reduction pro- 
cedure is to use the Kolmogorov backward equation to find a PDE corresponding to the 
system of SDEs we found. Specifically, given a sufficiently smooth function with compact 
support 6 : R N -> R, we define /(t,X,p,v) = £ x < p ' v [0(X t , p t , v t )]. Then 



at 



■■Mf), 



with initial conditions 



/(O,X,p,v) = 0(X,p,v) 



(7.1) 
(7-2) 
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Here, A e is the infinitesimal generator corresponding to the system of SDEs 4.10, 4.12, and 
4.13: 

Ae = L a l ow + -Lf ast (7-3) 



L slow ={-«-V(p - A[v]) ■ V x (p - A[v]) + [k- x C (p - A[v]) + A[-V X $(X)] (7.4) 

+ (V x A[v]) ■ v} ■ V p - C : V P + v ■ Vx 
Lfast =[-TC(v-v )] ■ V v +T : (7.5) 

Here, we use the double dot product, A : B = AijBij. We will use Einstein summation 
convention whenever a double index is present, unless otherwise noted. 

Next, we will perform a perturbative analysis on a part of A t . However, as we will see, 
to perform the perturbative analysis we will need a specific solvability condition to hold. In 
order to satisfy this condition, we will write 

A e = L 1 +L 1 + -L 2 . (7.6) 

and perform the analysis on the operator L £ = L\ + \L^ (we will take L2 = Lf ast ). Specif- 
ically, the solvability condition will require L\ to have mean with respect to the v distri- 
bution. The operator by which we will replace L\ + \L-z will be the first-order perturbation 
in e, while L\ represents the zeroth-order approximation in e. 

For situations with no slip (i.e. the limit e — > 0), we may use L\ as the effective operator 
replacing A t . Physically, we can think of v as being spread over a distribution in quasi- 
static equilibrium, and we may interpret L\ as the configuration average of L s i ow over the 
v distribution. That is, the fast variable v fluctuates so quickly that we may assume it is 
effectively in an equilibrium with p and X fixed. 

We can obtain the quasi-static equilibrium v-distribution, ip from the Kolmogorov for- 
ward equation, 

= [eLsiov, + £2]*V> (7-7) 

Taking e — > yields 

L* 2 ^ = (7.8) 

The solution to the above equation is a Gaussian, unique up to a multiplicative constant. One 
might expect to obtain such a distribution from Boltzmann-Gibbs statistics. Normalizing, 
we obtain 

^)=(^Pk e " |(V " V ° )TC(V " V0) - < 7 - 9 ) 

We will also use ip later in our perturbative analysis, although we will see analytically where 
the condition 7.8 arises. 

Taking the configuration average of v over L a i ow , we obtain 



Li = {-k V[P ' V x p - Av ■ V xP - p ■ V(Av ) + (Av ) • V(Av ) + (V X AC r ) : A T ] 

(7.10) 

+ K^Cip - Av ) + V X A : C- 1 + (VxA[v ]) ■ v - AV X $(X)} • V P + v • Vx - C : V P 

(7.11) 



The corresponding equations of motion are 
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^ = -k V[p ■ V x p - Av ■ V x p - p ■ V x (Av ) + (Av ) ■ V x (Avo) + (V X AC : A T ] 

(7.12) 

+ n^Cfa - Av ) + VxA : C- 1 + (VxA[v„]) ■ v - AVx$(X) + D-c^~ (7.13) 
^=v (7.14) 

Note that both V x and A are as a linear operators. This also allows us to compute the configu- 
ration average (using 7.27) (Av) • V(Av) = A i:j Vjd Xi A mn v n = d Xi A mn C~j Ajf = V X (AC _1 ) : 
A T . 

Here, we use V X A : C' 1 = [<9x fc A^C^ 1 ] and V X (A) 

re-labeling each variable with its physical dimensions (we assumed vo has dimensions of 
l/t 2 and that C has no dimensions) and using our original notation for gthm, we obtain the 
equations 

~ = -p~ V{p ■ V x p - A[mv ] ■ V x p - p ■ V x (A[mv ]) + (A[mv ]) ■ V x (A[mv ]) (7.15) 
at 

+ mfc s r(V x AC* _1 ) : A T } (7.16) 

+ P ~ X C (p - A[mv ]) + fc s TVxA : C' 1 + (V x A[mv„]) ■ v - AV X $(X) + A + g thm 

(7.17) 

§=v„ (7,8) 
Since we are assuming e — > 0, we may use vo = p _1 C _1 rp. 

7.2. e-Expansion of the Coefficients and Required Assumptions. Previ- 
ously, we found the PDEs that describe the evolution of the probability density function of 
the above SDEs using the Kolmogorov Backwards Equations. We will next isolate the op- 
erator that encapsulates the fast dynamics of the system that satisfies particular solvability 
conditions on which we will perform a perturbative analysis. 

In our reduction, we will replace the infinitesimal generator of a given stochastic process 
L E — L\ + by an effective first-order term in the limit e — > 0. We will require certain 
assumptions to hold for this reduction, which will be important in practice. 

To carry out the reduction, we consider the the stochastic process generated by ^ = 
L t f. First, we make a change of variable t = r/e to obtain 

®f = -L 1 + \L 2 (7.19) 
or e e 2 

and expand 

/(r, X, p, v) = /o (r, X, p, v) + e/i (r, X, p, v) + e 2 f 2 (r, X, p, v) + . . . 
We find the following relationships must hold by comparing powers of e: 

L 2 fo = (from \) (7.20) 

Lifa + La/i = (from -) (7.21) 

e 

^P- = Lift + L 2 f 2 (from 1) (7.22) 
or 

In our analysis, our choice of L 2 keeps X and p fixed. Assuming a system with L 2 as its 
infinitesimal generator is ergodic, we will now show that fo = fo(r, X, p). Since a solution 
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to 7.20 is a steady state of L2, we consider the PDE governing a function g(r, X, p, v) with 
any initial conditions g(0, X, p, v(0)) = #x, P (v(0)), 

L2 9 = IT- 
or 

In the limit r — > 00, g will reach a steady-state L^g = 0, which will be independent of v(0). 
By uniqueness, g ~ fo- 

We would like to use 7.21 to solve for /1. We notice a solvability condition from the re- 
quirement I/1/0 £ 7£(Z/2), or equivalently L1/0 G Af(L 2 ') ± ■ We can re-express this condition 
as follows: Vt/' £ N{L 2 ) (i.e. L^tp = 0), we must have 

(L 1 / Ol V)=0. (7.23) 

We find a 1/1 with the property j RN ipdw — 1 (given in 7.9). Checking the solvability condition 
with this choice of ip is sufficient since ip is unique up to a multiplicative constant. Supposing 
the solvability condition is satisfied, we can solve for 



From 7.22, we know 



Rearranging, we find that 



h = -L^Lxf . (7.24) 



-7^7 = —LiL 2 1 Lifo + Lifi- 



^- + L 1 L2 1 L 1 eN{Ll) 



Thus, 



(— +L 1 L 2 - 1 L 1 ,V) =0. 



But since fo is independent of v, 



/ 



dfo df f . , df 

-r—ipdv = — / ipdv = —. 



So 



dfo 
dr 



-[ 4)L 1 L 2 1 L 1 dv]f . 



(7.25) 



We use the resulting operator Lo = — [J & n ipLiL^ 1 Lidv] as the approximate generator 
to first order in e. After making back the change of variable from r to t, we see L ( is 
approximated by eLo in the limit e — > 0. Thus we replace 7.6 by the effective operator 

A = L 1 + eL (7.26) 

7.3. Preparing A e for the Reduction. Given tp, we must choose an appropriate 
L\ which will satisfy the solvability condition given by 7.23. Now, 



<L x f,i>>= \ L 1 (/)V = Li(/). 

Therefore, our choice of L\ must have configuration average 0. To achieve this, we split 
Lsiow appearing in 7.4 due to the slow variables as follows: 

Li L s iow L s i ow , L\ L s i ow 
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It will be useful to use the notation v = v — Vo. 
We will now compute V k Vj. 

First, we make a change of variable. We may assume that C, the covariance matrix, is 
positive-definite. Thus C has a square root matrix, which we may pick to be symmetric as 
well. Let v = C _3 ct. Noting the Jacobian ^f- = |C|~ 5 , we find: 



1 -1 -i _1 
(2n) N / 2 ° tj 2 C sk at0lae 2 1 l dai...da N 

C- h C; k hl = C^. (7.27) 



We then see 



V k Vj = (Vo k + Vk)(vOj + Vj) 



= VkVj + V k VQj + V 0k Vj + VQkVOj 
= C~k~j + VQkVQj. 

So that the configuration average 

(VxA[v]) • v = [d Xk A^v 3 v k ] = [d Xk KACk] + vo k v 0j )] (7.28) 
= VxA : C- 1 + (VxA[v ]) • v (7.29) 

From the above computation and v = we may obtain 7.10. 
to obtain Li, we compute 

VkVj - Wv] = v h Vj + v 0k Vj + VkVoj + v 0k v j - {v 0k v 0j + v k Vj) 

= VkVj + V 0k Vj + VkV j — VkVj 

It will be useful to maintain VkVj in L\ in its current form. 

For the purposes of a cleaner presentation, we shall use y = (X, p). 
We find Li is given by 

Q Q 

Li — Vi—— + [(— k~ CAv)i + dx k Aij(v k Vj + v 0k Vj + Vkv j — v k Vj)] — 

— k~ x v\— A„jVjd Xn Pi — p„dx n AijVj + A„ m v m d Xn AijV j + A nm v 0m d Xn AijVj (7.30) 
^ d 

+ A m kd Xm Aij(vkVj - v/cVj)]-^- (7-31) 
g 

= [RijVj + Sijk(vkVj — VkVj)] ——. (7.32) 
where 

R = ( -k-'CA + (V X A) ■ vo + Vx(Avo) - R 2 ) = ( B ) (7 ' 33) 

R 2 = «"V[-V x p ■ A - p ■ V X (A-) + V x (Av )A + Av • V X A] (7.34) 

S = ( V X A - k- 1 °[(V.A) ■ A T ] ) = ( P ) (7 - 35) 

The next step of the reduction, applying L^ 1 to Li, could be done term- by-term using 
elementary calculus. However, we prefer a more general and elegant approach, which we 
shall discuss next. 
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7.4. Inversion of L2. In summary, we will perform a change of basis on v so that 
L2 can be written as a sum of operators each involving only the i th coordinate. We will 
then show that L\ can be written as a linear combination of eigenfunctions of L2, each of 
which is a product of the eigenfunctions of the various L\ . 

7.4.1. Change of Variable. From 7.4 we have 

L 2 = -(TCv) ■ Vv + T : V| = -(TCv)i|i + T l3 -^ (7.36) 

OVi OViVj 

Since C is symmetric and positive-definite, we may find a square root C 3 which is also 
symmetric and positive-definite. By using a change of variable v = C~ 5 a, we may re- write 



We will use 

so that we may write 



L 2 = ClT l3 Cf m [~a m d an + d am d an ). (7.37) 

H = C5TC5, (7.38) 

L 2 = H nm [-a m 9 Qll + d am dc, n }. (7.39) 

Note H is symmetric since C 5 and T are symmetric, and 3 is positive-definite since T 
is positive-definite. Therefore, we may find a unitary decomposition S = QDQ -1 where D 
is diagonal and Q T = Q _1 . Note that in general we may have SJ 7^ A^C 3 . 

Now, take a = Qw to obtain 

L2 = D nm [~w m d Wn + d Wm d Wn ]. (7-40) 
To summarize, we have simplified the form of L2 using the change of basis 

v = C~3Q w = Qw. (7.41) 

7.4.2. Generic Formulation. We recall the Hermite equation — up'(w) +p"(w) = 
—Xp(w), whose eigenfunctions are the (probabilists') Hermite polynomials He\(w) with 
eigenvalues A corresponding to the degree of the polynomial. We see that for each i, the 
Hermite polynomials He\{wi) are eigenfunctions of (the term in L2 involving terms 
u>i) with eigenvalues —AD;. Thus TlkHe\ k (wk) is an eigenfunction of L 2 with eigenvalue 
— £\ XkDk- Moreover, any non-trivial product TlkHe\ k (wk) (i.e. at least some A*, 7^ 0), is 
in 1Z(L 2 ), since its corresponding eigenvalue is non-zero. 

Next, we will see that any polynomial p(w) that satisfies the solvability condition p(w) = 
can be written as a sum of non-trivial polynomials of the form YlkHe\ k (wk)- Consider 
the vector space Vi of such polynomials up to degree M. The solvability condition on p 
determines exactly one dimension, since a constant can be added to any polynomial to 
satisfy the solvability condition. 

Let V2 be the vector space generated by linear combinations of non-trivial polynomials 
of degree up to M of the form HkHe\ k Since the non-trivial Hermitian polynomials 

are orthogonal to constants, the non-triviality condition on the elements in V2 implies they 
satisfy the solvability condition WkHe\ k (wk) = 0. Thus V2 is a subspace of Vi. We see 
dimVi = Sd=i ( JV+ d _1 ) = dimV2, so we conclude Vi = V 2 - Thus we can express p(w) as 
a linear combination of non-trivial terms HkHe\ k (wk)- In practice, Gaussian elimination 
could be used to write p(w) in this form, starting with the highest-degree monomials and 
working downwards. 

Now, since a non-trivial term UkHe\ k (wk) is an eigenfunction of L 2 , we see 



L 2 1 [U k Hex k (w k ) 



YlkHe Xk (w k ). (7.42) 



Thus, noting L 2 1 is linear, the change of basis allows us to find L 2 1 p(w) 
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7.5. Computation for Perturbative Analysis. First, we re-write 7.9 in terms 
of w. Including the Jacobian, we obtain the distribution 

Next, we will compute 7.25. 

Plugging 7.41 into 7.32 and computing wjpw] using 7.43, 

Li = + S ljk (w k Wj - <5 fcj )J — = L^ 1 ' + Lf\ where (7.44) 

Rij — RikQkj , Sijk ~ SilmQljQmki (7-45) 

=R lj w j -^-, and if ) = SijkiwkWj - S kj )-^-. (7.46) 

Notice that each term in w is a Hermite polynomial, as we would anticipate. 

Some of the terms in our analysis will have zero contribution to the result. In particular, 
since i/>(w) is an even function in each of its coordinates, any polynomial in L\L~^L\ with odd 
multi-index will make a zero contribution to 7.25, since it must be odd in some component 
of w. Thus we may treat the terms L^ 1 and L^ 1 separtely. We will write 

li = - I ip^L^L^L^dw for i = 1, 2. (7.47) 

7.5.1. Computing I\. Using 7.42, integrating, and noting C (and therefore Q and 
D) are uniform in X, we compute 

h = - f ip^L^L^L^dw (7.48) 

tp(w)R nm w m dy n L2 1 [RijWjdy i ]dw (7.49) 

tp(-w)R nm w m dy n [RijDJ^Wkdy^dw = R n kD~^d Vn [Rijd yi ] (7.50) 

= RnrCrs 2 QskD^ 0y n [RipCpq 2 Qqj9 yi ] = Rnr^ r p ®Vn [Rip^Vi] (7-51) 
= Rnr^ rp 1 Rpidy nyi + T n p \dx n Bip] d Pi (7-52) 

= (Ry- 1 R T ) : Vy + [(VxB) : T" 1 ] ■ V p (7.53) 
= T" 1 : Vx + (BT' 1 B T ) : Vp + [(VxB) : T" 1 ] ■ V p (7.54) 

Above, we used that R is a function of X and not p, except in the first N rows where R is 
constant. 

7.5.2. Computing I 2 . Letting E mn = (D m +D„)~ 1 , using Sijk = S im ,„Q m jQ n k, 
and writing the sum over m and n explicitly to avoid confusion, 



h = - tp^L^L^L^dw (7.55) 

= — ^(w)Sijk(wjWk — Sj k )dy i L2 1 [Simn(wmW n — 5 mn )d yi ]dw (7.56) 
Jn 1 * 



In*' 

ip(w)Sijk(wjWk — Sj k )d yi 

R N 



dw (7.57) 



2, Sij k Si mn E mn / tj)(yv)(wjWk — Sjk) [(w m Wn — 5 )] dwd; m (7.58) 
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We note above Sijk only yields non-zero terms when i > N. 

An elegant way to evaluate the integral above (and more general similar integrals that 
may arise in a more general context) is to integrate by parts in one of the variables. Specif- 
ically, 



(2tt) jv/2 

1 



wje 2 Wj dwj Wk Wm w n ) e 2 ^ a ^ j Wa dwi...dwj...dwN 



(2tt) JV/2 

— & jk&mn ™f~ & jm&kn & jn&km 

where dvbj indicates a missing variable. 
Thus, 



_ 1 2 _ l 2 

e 2 ™i d Wj (wkW m Wn)dWje 2 ^ a ^j Wa dw\...dwj...dwN 



(7.59) 

(7.60) 

(7.61) 
(7.62) 



Finally, 



4>(w)(wjWk — 5j k )(w m w n — S mn )dw 

$ jk&mn ~\~ &jrn&kn H~~ $jn$km ' 2<5j/ c (5 m , n -|- Sjk^n 



E 

n^krn)dp- 

mn 

+ s 

d PiPi 



(7.63) 

(7.64) 
(7.65) 



(7.66) 
(7.67) 



Note for fixed m and n (no sum) we can also write 

mnSlnm — SiabSlcdQamQbni^QcmQdn + Q 

Using the form 7.68 in 7.67, we have 



^2 — SiabSlc 



E 

mn (QcmQdn + Q 



<~K, 



(7.68) 



(7.69) 



It may be useful to re-write 7.69 in the following way. We may interchange m and n. 
Averaging the original and new forms of I2 , we find 



I2 — ^SiabStcd 



+ Q anQbrn^Erri.n^QcrriQdn H~ QcnQd 



(7.70) 



This way, since the tensor applied to 6C. Pn is symmetric, we see I2 may be written as 



T 1 r T, a 2 It vv2 

h = ^\oa \ud p . Pl = -00 : V p 



(7.71) 



where it is a matrix. 
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7.6. Final Expression. We insert the results from the reduction into 7.26, using 
P = V X A to obtain 

A = Hf ■ V x + Hf ■ V p + Hf : V x + H 2 P : V P where (7.72) 

H 1 x =v , (7.73) 

Hf = -«rV[p • V xP - Av ■ V xP - p ■ V x (Av ( ,) + (Av ) ■ V x (Av ) + (V X AC _1 ) : A T ] 

(7.74) 

+ «r\C(p - Av ) + V X A : C- 1 + (VxA[v ]) ■ v - AV X $(X) + e(V x 5) : T" 1 

(7-75) 

Hf = eT" 1 (7.76) 
# 2 P = -C + e(BT- 1 B T + \o<J T ) (7.77) 

Next, we may write the equations of motion of the system (assuming the convention with p 
a continuum, but remaining in non-dimensional coordinates): 

* = Hf + S l hm (7.78) 

~=H? + Nl hm (7.79) 

where 

glbm = D H? £(X, t) and N^ m = D H xr)(t) (7.80) 
In the limit e ~ > 0, we obtain 

^ = -«"V[p ■ V xP - Ave • V x p - p ■ Vx(Av„) + (Av ) ■ Vx(Av„) + (VxAC -1 ) : A T ] 

(7-81) 

+ K^Cip - Av ) + V X A : C- 1 + (VxA[v ( ,]) ■ v - AV X $(X) (7.82) 
f=v (7.83) 

We may add back the units, understanding C has no dimensions and vo has units of i/t2, 
to obtain equations 5.1 and 5.2. 

8. Conclusions. An approach for obtaining systematic reductions of the fluid-structure 
equations in the limit of strong viscous coupling were developed. A dimension analysis of 
the stochastic fluid-structure interactions equations was performed. This analysis showed the 
importance of formulating the reduced equations in terms of the total momentum field of 
the fluid and immersed structures. The reductions were then developed based the Infitesmal 
Generators of the stochastic process and a singular perturbation analysis of the Backward 
Kolomogorov PDEs. Analysis was presented primarily in the physical regime where the 
viscous coupling that controls the momentum exchange between the fluid and structures 
became large. This established a few new SELM stochastic equations for modeling fluid- 
structure interactions taking into account inertial effects. We expect the SELM descriptions 
presented to be useful for a number of applications requiring accounting for inertial effects 
in fluid-structure interactions, such as in developing models for spectroscopy experiments. 
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